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Abstract 

We obtain Sobolev inequalities for the Shcrodinger operator —A — V, where V 
has critical behaviour V{x) = {{N — 2)/2)^|a;[^^ near the origin. We apply these 
inequalities to obtain pointwise estimates on the associated heat kernel, improving 
upon earlier results. 
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1 Introduction 

The purpose of this paper is to obtain some new Hardy-Sobolev inequalities and then use 
them in order to obtain new heat kernel estimates for the Schrodinger operator —A — V 
for positive potentials V with critical singularities, improving upon analogous estimates of 
this type. 

As a typical example, let us consider the case of a bounded domain f2 C R^, > 3, 
containing the origin. We obtain upper estimates on the heat kernel of the operator 

Hu = -Au - j^") ^1^!^ = 0) (1-1) 

for various values of the real parameter A. It is well known that the power is critical 
and that the associated heat kernel exhibits behaviour which is different from that of the 
case Ixl"''', 7 < 2, which is in the Kato class. Recently Milman and Semenov |MSj obtained 
upper heat kernel estimates for the operator when A < ((A^ — 2)/2)^. In the case 
< A < ((A^ — 2)/2)^ they showed that, modulo the usual Gaussian term, the heat kernel 
satisfies 

K{t,x,y) < ct"~|x|"°|?/|"". 
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where a denotes the smallest solution of aiN — 2 — a) = A. In Theorem 14.21 we extend 
this result to the critical case A = ((A^ — 2)/2)^, namely we prove that the corresponding 
kernel satisfies 

K{t,x,y) < ct-^'^\x\ —\y\ — —. 

This estimate is sharp as can be seen by comparing with the results of Vazquez and Zuazua 
|VZj . Moreover, for the subcritical case A < ((A^ — 2)/2)^ and for A < we improve 
the small-time dependence obtained in |MSj . namely we prove the estimate K(t,x,y) < 
ct~T+"|x|~°|?/|~" instead of K(t,x,y) < ct~~^'^~''\x\~°'\y\~°'; see Proposition 14.11 

In addition we consider operators that act on the whole of with potentials having 
the critical Hardy singularity near zero, of the form 



N-2 



2' 1-2 I I 1 

x\ , < 1, 



Veix) = { j i-^i ' p - (1.2) 

e/(x), \x\ > 1, 

under appropriate subcritical assumptions on the positive function /. Thus, in Theorem 
14. HI it is shown that if e > is small enough then the heat kernel of —A — K satisfies 

K{t,x,y) < ct"^/2max{|x|"^, l}max{|y|"^, 1}. (1.3) 

We also consider potentials that exhibit the critical behaviour ((A^ — 2)/2)^|x|~^ near 
infinity, that is 



e 



Under appropriate subcritical assumptions on g we obtain Sobolev estimates for —A — V^: 
for a sharp range of e > 0. We note here that while the question of Sobolev inequalities 
for —A — is rather similar to that for —A — V^, when it comes to heat kernel estimates 
essential differences arise. As mentioned earlier, the Sobolev inequality for —A — yields 
estimate ()1.3|1 for the corresponding heat kernel. On the other hand, while the short-time 
behaviour of the heat kernel of —A — is similar to that of the Laplacian, the long-time 
behaviour is very different. Zhang used a parabolic Harnack inequality to obtain precise 
estimates for the heat kernel of —A — V, when V is equal to A|a;|~^, A < ((A^ — 2)/2)^, 
near infinity. This generalizes the estimates given by Davies and Simon |DSj . where, 
however, more precise estimates were given. The corresponding problem for the critical 
case A = ((A^ — 2)/2)^ remains open. 

Going back to bounded Q C and to the operator H defined in fjl.l|) . for the critical 
case A = ((A^ — 2)/2)^ we finally consider additional singularities, that is, we consider 
potentials of the form ((A^ — 2)/2)^|x|^^ + Vi, where Vi > is also critical; Vi is defined 
as a series involving iterated logarithms (see definition (j5.17|) ) and is critical in the sense 
that the following improved Hardy inequality holds 

„|v.|^-(^)7^^>/^v^,.^ u.crm, (1.5) 

whereas this inequality is no longer true if we replace Vi by (1 + e)Vi for any e > 0. It 
is remarkable that the extra potential Vi does not affect the time dependence of the heat 
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kernel estimates, but only affects the spatial singularity at the origin (cf Theorems 14.21 and 
I5.3|) . This is in contrast with Proposition 14. If ii) where, for A < 0, the potential affects the 
the time singularity of the heat kernel as well. 

Throughout the paper we study a number of concrete potentials. These are chosen 
precisely because they are critical. By simple monotonicity one can then obtain heat kernel 
estimates for a whole range of other potentials, including potentials that are not radially 
symmetric. 

To prove the above heat kernel estimates we first use an appropriate change of variables, 
u = (pw, by means of which, the problem is reduced to obtaining uniform estimates on 
the heat kernel Kfp{t,x,y) of an auxiliary operator which acts on the function w; see, 
e.g., [MS]. Those estimates are in turn proved by means of some new Hardy- Sobolev 
inequalities. 

As a typical example of such an inequality we mention the following inequality proved 
by Brezis and Vazquez [BV]: 

valid for u E Hq{^1) and 1 < p < ]|32! this inequality fails for the critical Sobolev 
exponent p = -^Z2- To obtain sharp heat kernel estimates one needs to go up to the 
critical exponent. In connection with this we mention the following sharp Hardy-Sobolev 
inequality established in [FT]: 

valid for u G HliVL); here D = sup^ |x| and Xi{t) = (1 — logt)~^, t G (0, 1). In the present 
work we derive new Hardy-Sobolev inequalities that involve potentials such as the ones 
given in ()1.2|) or ()1.4|) : see Theorems 13. 4| 13. 5| 15.11 and 15.21 We should mention that the 
validity of improved Hardy inequalities is strongly connected to the existence and large 
time behaviour of solutions of the heat equation with singular potential; see, e.g., Brezis 
and Vazquez [B V] , Cabre and Martel [CM] , Davila and Dupaigne [DD] as well as Vazquez 
and Zuazua [VZ]. 

As a byproduct of our approach we establish various results concerning improved Hardy 
inequalities with boundary terms. Such inequalities have recently attracted attention, see 
Adimurthi and Esteban [AE], Wang and Zhu [WZ] and references therein. 

The structure of the paper is as follows: in Section 2 we present some auxiliary results 
concerning improved Hardy inequalities with boundary terms. In Section 3 we prove the 
Hardy-Sobolev inequalities; in Section 4 we apply them to obtain heat kernel estimates. 
Finally, in Section 5 we prove refined Sobolev inequalities and heat kernel estimates when 
additional singularities are present. 

Acknowledgment We acknowledge partial support by the RTN European network 
Fronts-Singularities, HPRN-CT-2002-00274. 
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2 Two minimization problems 



Throughout this section, C R , > 3, is a bounded domain containing the origin with 
boundary. Also, we always denote by u the outward-pointing (with respect to Q) unit 
vector on the surface dfl. In Section 2.1 we will work on Q, while in Section 2.2 we will 
work on fl"^. The results of this section will be applied in Section 3. 



2.1 Bounded domains 

For a > we define 

Ma) = inf — ^ ^^^^ . (2.1) 

Lemma 2.1 Let Q be a bounded domain in R^, iV > 3, containing the origin. 

(i) If < a < then Xn{a) = a{N — 2 — a). Moreover, \x\~°' G H^{Q) is a minimizer 
for < a < ^^^Y^, whereas for a = there is no H^{^1) minimizer. 

(ii) If a > ^^-^ and Vt is starshaped with respect to zero, then Xn{a) = ("^^^y^) o-i^d th 
is no II^{Q) minimizer. 



ere 



In case Q is not starshaped with respect to zero, concerning the analogue of part (ii) 
of the above Lemma, we have 

Lemma 2.2 Let Q be a bounded domain in R^, N > 3, containing the origin, which is 
not starshaped with respect to zero. Then, there exist finite constants a* > N — 2 and 
a* G [^^^,a*) depending on f2 such that: 

(i) \n{a*) = 0, whereas Xn{a) > for all < a < a* . 

(ii) // ^^Y^ < a < tt*, then \q,{q) = {^^-^) and there is no H^iVt) minimizer. 

(iii) If < a < a*, then max(0, a{N — 2 — a)) < Xn{<y) < i^^^)'^ , and there exists an 
H^{fl) minimizer. 

Remarks. 1. We note in particular that for any Q and any a > there holds 

a{N-2-a)< \n{a) < (^— ) • (2-2) 

2. We do not know whether there exists a non-starshaped domain O with smooth boundary 
so that a* = N — 2. Similarly, we do not know whether there there exists such an Q for 
which a, = {N -2)/2. 

Proof of Lemmas \2.1\ and \2.i^ Let u G C°°(i7) be supported outside a neighborhood 
of zero. For any a > we set u{x) = (x). A straightforward calculation shows that 



/ \Vu\'^dx= I |x|-2°|Vfprfx + a(A^-2-a) f ^dx - a f ^^u^dS, (2.3) 
Jn Jn Jn \x\^ Jan |a;p 

therefore, 



\Wuydx + a --:rU^dS>a(N-2-a) —dx. (2.4) 
JdQ \xr Jn \xr 
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By a simple density argument this inequality is valid for all u G H^{Q). This implies in 
particular the lower bound on An (a) in ()2.2p . 

If < a < then | x\ " is in H^{il) and an easy calculation shows that satisfies 

()2.4|1 as equality, hence it is a minimizer. If a = then the fact that \n{a) = ("^^2^) 

follows by considering the minimizing sequence u^{x) = \x\ 2 '"^ G H^{Q), e 0^. 

By the same sequence Uf{x) = Ixl""^"^^ G H^{^1) one can show that ("^^f^) > An(a), 
for any a > and any Q, thus proving the upper bound in ()2.2j) . 

Suppose now that f2 is starshaped and a > Then, using first the fact that x-u >0 
on the boundary of Q, and then ()2.4|) (with 



a 



2 



Wu\''dx + a I ^^u^dS > f iVul^dx + ^^—^ f %^u^dS 



an \xr 2 Jan \x 

N-2^2 f 



> { — ) / r^dx. 

Jn Ixr 



Hence, in this case Xn{<y) = {^^y^Y ■ 

We next show that when XQ,{a) = (^^^y^)^, there is no H^{Q) minimizer. Indeed 
assuming that there is one, then it would be a positive H^{Q) solution of the Euler- 
Lagrange equation 

An+ u = 0, xen, 

\x\'^ 

under suitable Robin boundary conditions. However, this equation has no H^{Q) positive 
solutions; see e.g. [FT, Theorem C] for a more general statement. (We note that although 
in this Theorem Dirichlet condition were imposed, the proof is independent of the boundary 
conditions.) Thus, Lemma f2. II has been proved. 

Suppose now that Q is not starshaped with respect to zero. The existence of a* follows 
from the continuity of Xn{a) with respect to a combined with the fact that if fl is not 
starshaped with respect to zero, then one can easily find test functions making the surface 
integral in ()2.ip negative. The fact that a* > N — 2 follows from the lower bound in ()2.2p . 

From Lemma 12.1( 1). we have that X^{^) = {^f. We then define a^: as the 
supremum of all a for which Xn{a) = i^^)'^- Assuming that a* > we will show that 
for any < a < a* there holds Xn{a) = {^^y^Y ■ Indeed, if this is not the case then 
there would exist an a in the above interval and cj) G H^{Vt) such that 



< (2-5) 



On the other hand from Lemma [2.1( 1). we have that 

\VcP\'dx + (^) J,^ f^^'dS ^ /AT - 2^ 2 

Inj^dx - V 2 ^ ^ - > 

From the above two inequalities it follows that Jqq ^(p'^dS < 0. Using as a test function 

and the fact that a^, > a we conclude that X^{a^) < {^^^f, which is a contradiction. 
Thus, the estimates of part (ii) and (iii) of Lemma 12.21 have been proved. 
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The nonexistence of H^{Q) minimizer of part (ii) follows exactly as in Lemma ITTl The 
existence of H^{Q) minimizer of part (iii) will follow later from a more general result; see 
Proposition 12 .51 / / 



2.2 Complement of bounded domains 

Here we consider the complement of a bounded domain and we study the corresponding 
infimum, that is 

Un(a) = mf 5 • (2-7) 

where fi, as before, is a bounded domain containing the origin and v is the outward- 
pointing (with respect to Vt) unit vector on the surface dVL. Also, C^(R^)|n= is the set of 
restrictions on Vf^ of all functions u G C^(R^). We also introduce the following norms, 

||M||^i,2mc) = (/ |VMprfx)^/2 + (/ Iwl^rfx)^ (2.8) 



if \Vu\^dxy/^ + { f ^rfx)^/2 (2.9) 



/ o 1 /n I 2(JV-1) N-2 

u\\yv(n^) = (/ |VMpdx)^/2 + ( / \u\^^dS)^^ (2.10) 



dn 

and we denote by ^'^'^(fi^), n\Q^) and W{Q^) the completion of C^{R^)\nc under the 
corresponding norms. The space W(fi^) is well studied, see e.g. IMj- For our purposes 
however, the natural spaces to use are P^'^(f2'^) and T-C^{Q'^). In the next lemma we show 
that these three spaces coincide (a trivial fact if Q'^ were replaced by Q). 

Lemma 2.3 Let Q be a bounded domain in R^, > 3, with boundary, containing the 
ongm. Then ©^'^(fi^) = n^iVt^) = W(fi'=). 

Proof: We will show that all norms are equivalent. Let u G C^(R^)|qc. Under our 
assumptions, it follows easily as in Lemma 0(cf (Q with a = ^) that 

^u\^dx — / -—u^dS>{^—f —dx, (2.11 

ri'^ ' 2 Jan \xr 2 Jn'^ \xr 



whence. 



-dx < C 



Q= \x\ 



( I \Vu\^dx+ I \u\'^dS 
^Jq,'^ Jan 



< cU \Vu\^dx+U \u\^^dS 

1/ 1/ 



N-2 
(JV-1) 



/an 

Hence, ||^i||-Hi(nc) ^ C'||^||w(n=)- To obtain the reverse inequality we note that it follows 
from the standard trace Theorem (e.g [A, Theorem 5.22], Chapter V) - applied to B \ Q 
for some ball B D Q - that 



(Jjul'^dS)^ <C\\u\\nHn^^ 



u\'-^dS < Cllnlki 



an 



From the first one it follows that ||'u||w;(r2c) < C||m||7^i(qc), whence H^{fl'^) = W{Q'^). From 
the second one it follows that ||M||w(r2<:) < C||w||x)1,2(qc). Thus, it remains to prove that 
||m||x)1,2(qc) < C||m||vv;(qc). This inequality follows from Corollary 1 of Section 4.11.1 [M, p. 
258]. Notice that in the notation of Maz'ja W(fi^) = iy„ 2(iv-i) ((^^ dn). // 

iV-2 

An immediate consequence of the above Lemma is that the infimum in ()2.7p can be 
taken over T>^''^{Q'^), that is 

Sac \Vu\'^dx-a Jga rriU^dS 
Mo^)= inf T^^^^ . (2.12) 

We now state the analogue of Lemma f2. II for exterior domains. 

Lemma 2.4 Let Q be a bounded domain in R^, > 3, containing the origin. 

(i) If a > ^^Y^, then fiQ^a) = a{N — 2 — a). Moreover, G is a minimizer 
for a > ^^Y^, whereas for a = there is no minimizer. 

(ii) I/O < a < ^^^Y^, and f2 starshaped with respect to zero, then fin{a) = ("^^y^) there 
is no V^''^{Q'^) minimizer. 

Proof: The proof is quite similar to the proof of the previous Lemmas 12.11 and 12.21 An 
alternative proof can be given using the Kelvin transform; see the Remark that follows. / / 

Remark There is a duality between the minimization problems ()2.1|) and ()2.7|) . Indeed, 
by means of the Kelvin transform, u{x) = \y\'^^'^v{y), y = x G Q'^ the domain Q'^ is 

transformed to a bounded domain containing the origin that we denote by [Q'^)*. Denoting 
by u* the outward pointing normal to d{Q^)* a straightforward calculation shows that 

/ \VM^dx= f \Vyv\^dy+(N -2) f ^^-^v^dSy. 

Jn-^ Jim- Jdiu'^)* y r 



Also, 



T^dx = / Y^dy, 

2JV f 2N 

lul^-'^dx = / \v\'^~^dy. 



It can be seen from these relations that u G T>^''^{Q'^) if and only if f G H^{{Q^)*). It then 
follows easily that /in(«) = A(oc)*(A^ — 2 — a), and that the existence of a minimizer for 
/in (a) in is equivalent to the existence of a minimizer in H^(fi) for A(nc)*(A^— 2— a). 



2.3 Existence of minimizers 

In this section we establish a sufficient condition for the existence of minimizers. We recall 
from Lemma f2. 21 that when Q is not starshaped with respect to the origin, a* denotes the 
first zero of An (a). We also set a* = oo in case Q is starshaped with respect to zero. Thus, 
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in both cases we have Xn{<y) > for < a < a*. Given < a < a*, and a nonnegative 
measurable potential V we define 

Mo..V):= inf , . (2.13) 

Note that with this notation An(«) = Aq(q;, |x|~^). Since the numerator in ()2.13p is always 
positive and finite when < a < a*, we interpret \^{a,V) = in case there exists 
u e H^iyt) such that j^Vv?dx = +oo. It is worth mentioning that Xfi{a,V) is not 
monotone with respect to Q, unlike the case of Dirichlet boundary conditions. 

We denote by Bj. C Q the ball centered at zero with radius r. We have the following 

Proposition 2.5 Let Q be a bounded domain in R^, > 3, containing the origin, and 

JV 

letO<V e Lil^{n \ {0}). If for some r > 

< Xn{a,V) < XBAa,V) (2.14) 
then has an H^{VL) mimmizer. 

Note. It is a consequence of ()2.14|1 that Vv^ < +oo for u G H^iVt). 

Proof. Let {uj} G H^{VL) be a minimizing sequence of the Rayleigh quotient in ()2.13j) . 
We may normalize it so that Vu'jdx = 1. We claim that < C. This will follow 

from two inequalities. The first inequality follows from the fact that < a < a* and 
Xn{(^*) = and reads 

Vu\^dx + a[ ^u^dS > {1 - —) f \Vu\^dx. (2.15) 
n ' Jan \xr a* Jn 



The second one is a consequence of Lemma 12.11 and reads 



Vu\ dx + a , ,^ u dS > X^ia) / ■; — —dx 



n Jan \xr Jn x'^ 



> KXn{a) / u'dx. (2.16) 
Jn 

Thus, we may extract a subsequence such that Uj Uq weakly in H^{VL), and Uj — >■ Uq 
strongly in -L^(f2), 1 < p < Moreover, since V G L^^'^{VL \ Br), standard results give 



Also by the trace theorem 



Vu^.dx-^ / Vuldx. (2.17) 

n\Br Jn\Br 



^u]dS^ I "^uldS (2.18) 
an \xr Jan \xr 



Setting Uj = vj + uq we easily see that as j — oo, 

J^\\/Uj\'^dx = J^\\/uo\'^dx + J \Vvj\'^dx + o{l), (2.19) 



and 



1= / Vu'^.dx= f Vuldx+ f VvUx + o(l). (2.20) 
Jn Jn Jn 

It then follows from ()2.13p that 

Xn((y,V) = / iVvjl'^dx + / I'Vuofdx + a / —rruldS + o(l) 
Jn Jn Jan |xp 

> / iVvA'^dx + Xn(a,V) [ Vuldx + o(l). (2.21) 
Jn Jn 



We then have 



|Vfj| dx > 



Br 





dx 


a,V) 




a,V) 




a,V) 


Jn 



vIdS 



= XBAa,V)[l- J^Vuldx)+o{l) (j^oo). (2.22) 

Using this and ()2.2H) we end up with 

(An(a, V) - \B,Xa, V)) (l - ^ Vuldx) > 0, (2.23) 

whence, since Ao(a, V) < Xs^^a, V), it follows that Vu^dx > 1. By lower semi continuity 
we conclude that Vu^dx = 1. It then follows that Uq is a minimizer for ()2.13|) . / / 

As a consequence we have: 
Completion of Proof of Lemma \2.^^ iii) (Existence of a minimizer) : Since ^^-^ < a* < 
a < a* it follows from Lemma O that < X^{a) < {^f. liBrCVt is a ball centered at 
zero it follows by Lemma ITT] that AB^(a) = {^^-^)^ . By Proposition 12.51 An(«) is attained 
by an H^iVl) function. // 

We next state the corresponding result for the exterior of a bounded domain Vl. For 
0<a<A^-2we define 

^in{a,V)■.= inf ' , , I"' . (2.24) 

^ ^ ' ' uev^-Hn-) J^.Vu^dx ^ ' 

> 

We then have 

Proposition 2.6 Let Q be a bounded domain in R^, > 3, containing the origin, and 

N_ 

letO <V e LiI^{VL''). Also let < a < N - 2. If for some ball BjiD centered at zero 

0</if,(a,\/) <^B«(«,y), (2.25) 
then \2. 24\ ) has a minimizer. 
The proof is similar to that of the previous Proposition. 
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3 Hardy-Sobolev inequalities 
3.1 Auxilliary inequalities 

We begin this section with two known Sobolev inequahties that will be used in the sequel. 
In Theorems IH.4I and IH.5I we then prove two new Sobolev inequalities. 

By the classical inequality of Caffarelli-Kohn-Nirenberg [("KNj we have that 

J^^\Vw\^\x\^^''dx>c[J^^\w\P\x\-P'^dxf^\ weC^{R^), (3.1) 

with p = 2N/{N -2 + 2{(3- a)), provided a < (A^ - 2)/2 and < /3 - a < 1. 

At the critical case a = P = {N — 2)/2 inequality ()3.1|1 fails. A sharp substitute for 
bounded Q was obtained in |FTj . where it was shown that, with 



Xi(t) = (l-logt)-\ tG(0,l), 
and D = sup^ \x\ there holds 

/ \\/w\'\x\'-''dx>c([ (3.2) 

for all w G C^(f2), where the exponent of Xi(|a;|/D) is optimal. 

In the sequel we will make essential use of the following one-dimensional result, which 
is a special case of a more general statement by Maz'ja, cf. p31 Theorem 3, p. 44]: 

Proposition 3.1 Let Air), B{r) nonnegative functions such that 1/A{r) and B{r) are 
integrable in {r,oo) and (0,r) respectively, for all positive r < oo. Then, for q > 2 the 
Sobolev inequality 

r°° / \ 2/a 

{v'{r)fA{r)dr>c[j^ |t;(r)|''5(r)rfr) , 

is valid for all v G C^(0, oo) that vanish near infinity, if and only if 

dt ' 

r>0 ^"'0 ' A{t)' 



sup { / B{t)dt) { / -^Y'"^ < 



The above proposition will be applied to higher dimensions by means of the following 

Lemma 3.2 Let N > 2. Suppose that V G ^^^^(R^ \ {0}) n LL(R^) «s a radially 
symmetric function. We further assume that inequality 

[ IVmPc/x - / Vu^dx > 0, (3.3) 
is valid for all radially symmetric functions u G C^(R^). 

(i) Then, \3. is also valid for nonradial functions, that is, for all u G C^(R^) 

(ii) //, in addition, 

< ess sup^gj^jv |xpV^(x) = 6 < GO, (3.4) 



10 



then the following improved inequality holds 

iVuPcix - / Vu'^dx > (3.5) 

Jrn 



> 



f \Vuo\^dx- [ Vuldx + , \ / \V{u~Uo)\'^dx 



where Uo{r) denotes the spherical average of u & C^(R^), that is 

NuJNr^~'^ JdBr{0) 

Proof Let u E C^(R^) and let 



^o{r) = I u{x)dS^. (3.6) 



oo 

'2f cN~l\ 



be its decomposition into spherical harmonics; here fm are orthogonal in L (S* ), nor- 
malized by /5JV-1 fi{a)fj{a)dS = 6ij. In particular /o(cr) = 1 and the first term in the 
above decomposition is given by (j3.6|) . The /^'s are eigenfunctions of the Laplace-Beltrami 
operator, with corresponding eigenvalues Cm = m{N — 2 + m), m > 0. An easy calculation 
shows that 

I.dV-P - Vu^)dx = E /^J|V-.P + (^ - V)ul}dx = 

= / (|V«oP - Vul)dx + / J|V«„.P + (fi - ^)«^}^^- (3-7) 

Part (i) follows immediately since > and radially symmetric. 

To prove part (ii) we first observe that 

j^^ I V(n - Uo)\Mx = J2j^^{\^^rn\^ + J^^^l}d^- 
In view of this and ()3.7|1 it is enough to establish that for any m > 1, there holds 

or, equivalently, 

/ \Vum\'^dx> [ ul^\— — ^^-^V - -^}dx. 
Jrn ' Jrn ™L \x\^) 



Since > Ci = — 1 it is enough to establish this for Cm = N ~ 1. By the definition of 
6, cf fl3.4j) . it follows easily that 

and the result follows from (j3.3p . / / 

As a consequence of this we next establish the following result. 
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Lemma 3.3 Let N > 3. Suppose that V e Iv^,(R^ \ {0}) n Ll^iR^) ts a radially 
symmetric function, such that 

< ess sup^gpj^jv = < oo, 

and W G L°°(R^) is a positive radially symmetric function. We further assume that the 
inequality 

/ \Wu\^dx- Vu^dx>c( (3.9) 

is valid for all radially symmetric functions u G C^(R^). Then inequality 4,'/. .9)) is true 
for all u G C7^(R^) (without radial symmetry), provided the constant c is replaced by a 
new constant C depending on c, N, 6 and \\W\\loo. 



Proof: Starting from ()3.5j) we compute 

J^^ {|Vm|2 - Vu'^}dx > 

f 1\ — 1 f 

/ r \ iN-2)/N 

> ci^J^Juor/^''-'^Wdxj + 

> C( f 



N-1 
-1 + 

{N-2)/N 

{N-2)/N 
{N~2)/N 



where, for the last inequahties we used the standard Sobolev inequahty, the boundedness 
of W and the triangle inequality. / / 

3.2 Hardy-Sobolev inequalities 

In this section we prove improved Hardy-Sobolev inequalities for potentials that are critical 
either near zero or near infinity. We first consider a potential which is critical near zero. 
For e > we define 



K(x) = { V 2 J i-i ' - - ^ ^ (3.10) 
e/W, \x\ > 1 

where / is a non-negative, continuous and radially symmetric function on {|x| > 1}. 
Moreover we assume / to be subcritical, satisfying 

f{x) < K\x\~^~'', \x\>l, (3.11) 

for some a,K > 0. 

We also define the following auxiliary function (cf. ()5.15|) ) 



Xi(|x|)= f^^l-^l^' ^ ■ (3-12) 





\x\ 


1 ^l(l^l), 


\x\ 
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We shall henceforth denote by B the unit ball in centered at zero, by i?^ its 
complement, and, as before, we denote by C^(R^)|bc the set of restrictions on of all 
functions u G C^(R^). We also denote by u the outward-pointing unit vector on the 
surface dB. We have the following 

Theorem 3.4 Let 

Then eo > and for any e G (0, eo) there holds 

j^^ \Vu\^dx - J^^ V^u'dx > c (/^^ \ur'^''-^^xf''-'^'^''-'\\x\)dxY~'^''' , (3.14) 
for all u e C^{n^). Moreover, fails for e = eo- 

Proof. Since f{x) < K\x\~^ the positivity of eo follows from Lemma 12.41 with a = 
yielding in fact eo > K-\{N - 2)/2f. 

Let us now fix e G (0, eo). By Lemma ITS! it is enough to prove ()3.14|) in the case where 
u is radially symmetric, u = u{r). Now, there exists a radially symmetric and positive 
function if) on B'^ which solves the Robin problem 



A?/^ + e/V' = 0, |x| > 1 

dv 



The existence of such a ip can be easily derived, for example, by a shooting argument from 
= 1}. We assume that ijj is normalized so that ip = 1 on {\x\ = 1} . The function 




then lies in C^(R^\{0}), is positive, radially symmetric and satisfies Aip + Veip = in R'^. 
Following JFT' we change variables, u = ipv, and ()3.14j) for radially symmetric functions 
is then written as 

/ rOD \ IN -2) IN 

^ (t;') Vr^^'^r > C |^|27V/(^-2)^27V/(^-2)^2(Ar-l)/(Ar-2)^^y _ ^3_^g^ 

We claim that %l){r) has a positive limit as r — > +oo. Indeed, since {r^~^tp')' = — r^~^ V^?/' < 
and ip'll) < 0, ilj{r) is decreasing on (l,+oo). If the limit lim^^+oo were zero it 
would then follow from |LJN| Theorem 2.9] that ?/'(r) < cr^~^ near infinity, which then 
easily implies G 1H}{B'^). Hence ip can be taken as a test function for the infimum in 
the right-hand side of ()3.13p . in which case the value of the Rayleigh quotient is e < eo, 
contradicting the definition of eo. Hence lim^^+oo V^(^) = / > 0. Using this we deduce 
()3.16|1 from Proposition 13.11 Hence ()3.14|1 has been proved. 

We finally show that ()3.14j) fails for e = eo. For this we will use Proposition 12.61 Let 
BrD Bi. Then V,,{x) < eoKR-''\x\-^ on {BrY, hence 
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By Lemma EiH (i), for a G (0, — 2) /iB^(Q;), is positive and independent of R; taking 
R large enough we have /iB(tt,V^o) < /^Bfllc^^^o)- Hence, by Proposition 12.61 - with 
a = {N — 2)/2 - there exists an ?-^^(i?'^)-minimizer to ()3.13p . It is standard to show that 
is simple, radial and of one sign; we normalize it by 0||a;|=i = 1 and for > we define 
the function ug G H^(Il^) by 



ue[x 



I I ^-2 I I I 

\x\ 2 \x\ < 1: 



[ 0(a;), |x| > 1. 

We then compute the left-hand side of ()3.14|) : in B there holds Auq + V^^gUg = hence 

{\Vug\^ -V,yg)dx = 



- / {ugAug + VeoUl)dx + [ Ug-^dS + [ {\Vug\^ - V^^uDdx 
Jb JdB ov Jb^ 



B r"^ '2 '2 JdB 



_ Nuon9 
2 ■ 

On the other hand for the right-hand side of ()3.14|) we have 

the last term being independent of 6. Letting ^ we conclude that ()3.14|) fails for 

e = eo- // 
We close this section proving a Sobolev inequality which involves radial potentials 
with critical behaviour ((A^ — 2)/2)^|x|~^ near infinity. Let g he a non-negative, radially 
symmetric, and continuous in i? \ {0} function that is subcritical near zero, satisfying 

g{x) < K\x\~^+'', \x\<l, 

for some a, K > 0. For e > we define 



We also set 

We then have 
Theorem 3.5 Let 



Yi{\x\) 



L, |x| < 1, 

1 + In \x\y^, \x\ > 1. 



. J^\Wu\'dx + ^Js^u'dS 
eo = mf 



ueH^(B) J^gu'^dx 
Then eo > and for any e G (0, eo) there holds 



I \Vu\^dx - I V^u'dx >c( f |np^/(^-^)y;'(^-^)/(^-')(|x|)da;V''~'^^'^ , (3.18) 
Jr'^ Jr^ \Jr^ \i 1/ y 

for all u G C'^{BJ^). Moreover, HH^) fails for e = eo- 
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Proof. The proof of ()3.18|) follows closely that of Theorem 13 .41 reversing essentially the role 
of B and B'^ while making the necessary adjustments; in particular, we now use Lemma l^?T] 
instead of Lemma (2.41 In fact, an alternative and simpler proof consists in simply taking 
the Kelvin transform of (j3.14|) . The optimality of eo is also proven analogously; we omit 
the details. / / 



4 Heat kernel estimates 

In this section we shall apply the Sobolev inequalities of Section 3 to obtain heat kernel 
estimates for the Schrodinger operator 

Hu = —Au — Vu, u\dn = 0, 

for various critical potentials V. We always assume that is a domain containing the origin 
in R^, > 3. We shall consider the case of bounded Q, as well as the case Q = R^. The 
operator H is defined via quadratic forms, with initial domain C^(fi \ {0}); it will always 
be the case that H > 0. Note that, equivalently, we could have set Cl{Q) as the initial 
domain. 

We shall use the standard technique of transference to a weighted space, which we 
now describe briefly. Let (j) E C^{Q\ {0}) be positive and such that A0 G Lj^ci^ \ {0})- 
The unitary map 

L\n)3u^w = ^eLl:=L\n,(l)^dx) (4.1) 



satisfies 

/ (I Vwp - Vu^)dx = f ( |Vwp - Vw^ - ^w^)(t)Hx (4.2) 
Jo. Jn^ (p 

for all u & CI{VL \ {0}). Hence, if in addition satisfies A0 + Vcf) = (weakly) on i7 \ {0}, 
then 

/ (iVwp - Vu^)dx = I \Vw\^(t)'^dx (4.3) 
Jn Jn 

for all u G C^(f2\{0}). Hence H is unitarily equivalent via ()4.1|) to the self-adjoint operator 
on L^, defined initially on C^(fi \ {0}) and given formally by 

H^w = -^div(0^VM;), = 0. 

The space Cl{Q \ {0}) is invariant under multiplication (or division) by and hence it is 
a form core also for if^. Moreover, a Sobolev inequality of the form 



/ (|Vm|2 - Vu^)dx >c(f \u\''Wdx 
Jn \Jn 

is valid for all m G C^(r2 \ {0}) if and only if 



2/q 



2/q 

iVwiVc/x > cf / \w\''4>''Wdx^ 
n 
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for all w ^ Cl{Q\ {0}). Finally, the heat kernels of H and are related by 

K{t, X, y) = (t){x)(l){y)K^{t, x,y), t>0, x,yen, (4.4) 

and hence one can obtain estimates on K{t, x, y) via estimates on K^{t, x, y). 

Example. As a typical example let us consider the case of a bounded domain Q in R^, 
N >3, and let V{x) = \\x\-'^, \ < {{N - 2)/2)2. Let (j){x) = a being the smallest 

solution of a{N — 2 — a) = X. Then A0 + V(j) = on Q \ {0} and therefore 

/ (|VMp-A^)rfx= / |Vwnx|-2°cix. (4.5) 
Jn |a;p Jn 

for all u ^ Cl{Q \ {0}) or, equivalently, for all w & Cl{Q \ {0}). Moreover, the heat kernel 
of H = —A — V is related to the heat kernel of by 

K{t,x,y) = \x\-°'\y\~°'K^{t,x,y). 

We note here that a simple approximation argument shows that for A < ((A^ — 2)/2)^ the 
form domain of H is Hq{Q), but at the critical case A = {{N — 2)/2)^ the form domain is 
strictly larger than H^{fl); see also |FTj . 

Sobolev inequalities are related to heat kernel estimates by the following standard 
result jDj Theorem 2.4.2]: for any q > 2, 



the upper bound 

K^{t, X, y) < ct"*?/^, t > 0, x,y e Q 
* is equivalent to the Sobolev inequality 

(4.6) 

J^\Vw\'<f>'dx>c{j^\w\^-<j>'dxy'''^^\ weCl{n\{0}). 

In the rest of this section we shall apply the Hardy-Sobolev inequalities of Section 3 in 
order to obtain upper estimates on the heat kernel K{t, x, y) of the operator —A — V for 
critical and subcritical potentials V. For this we shall use (j4.4j) for appropriate functions 
(f), together with uniform estimates on K^(t, x,y), obtained by means of (j4.6|) . We present 
on-diagonal estimates, in the sense that we do not include the usual Gaussian term, which 
can be added by standard methods. 

We assume that f2 is a domain in R^, > 3. We retain the notation introduced 
in the last example, and, in particular, we have H = —A — A|x|~^, subject to Dirichlet 
boundary conditions on dQ. Part (i) of the next proposition has been proved in |MSj but 
we include a proof for the sake of completeness. Part (ii) improves upon the corresponding 
estimate of |MSj . who had K{t,x,y) < ct 2 ''la;!""!!/!"". 

Proposition 4.1 (subcritical case) Let K(t, x, y) he the heat kernel of H = — A — Xj^2, 
subject to Dirichlet boundary conditions on dQ. For A < ((A^ — 2)/2)^, let a be the smallest 
solution of a{N — 2 — a) = X. 
(i) Ifn IS bounded and < X < {{N - 2)/2f then 

K{t,x,y) < ct~^/'^\x\-''\y\-'', t>0,x,yen. (4.7) 
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(ii) For any Q C (bounded or unbounded) and A < then 

Kit, X, y) < ct-""^ mm{l, (^)-n min{l, (^)-"}, t > 0, x, y G fi. (4.8) 

Proof, (i) The boundedness of together with imply 

^ iVwHxr^^rfa; > \w\-^^\x\-^''dx)''^'^^'^ , w G C^{n). (4.9) 

By ()4.fjj) this imphes K(p{t,x,y) < ct~^^'^, from which ()4.7|1 follows using ()4.4|1 . 

(ii) Comparison with the Laplacian implies that K{t,x,y) < ct~^^'^. Moreover, in- 
equality p.ip for Pp = 2a reads 



/ 

Jn 



2(N-2 ) \ -'^-2-2 



By means of (|4.6p we deduce that K^{t,x,y) < ct~^/^+°, t > 0, x,y E ^l. Hence 

K{t,x,y) <ct-''/'mm{l,{^)-^{^r^}, t>0,x,yeQ. (4.10) 
This proves ()4.8|) when x = y. The general case follows from the semigroup property since 



K{t,x,y) = I^K{t/2,x,z)K{t/2,z,y)dz 

1/2/ r ^ , ^„ \ 1/2 



n 



< 



K{t/2,x,zfdz) (/ K{t/2,z,yfdz) 
ir(t,a;,x)i/2ir(t,l/,y)V2. 



// 



Theorem 4.2 (critical case) Lei Q be a bounded domain and K{t, x, y) be the heat kernel 



of H = —A — ("^^2^)^^; subject to Dirichlet boundary conditions on d^l. Th 



en 



N , , JV-2 , , JV-2 



K{t,x,y)<ct 2 |x| 2 \y\ 2 ^ t>0,x,ye^. (4-11) 

Proof. Estimate ()3.2j) implies the weaker inequality 

1^ iVwl'lxl'-'^dx > c[j^ It/^li^lxl^-^rfx)^"^"'^^"^. (4.12) 



Hence K^{t,x,y) < ct~^ and ()4.11|) follows. // 
We next consider the case Q = R^, and the potential is critical at zero. More precisely 
we consider the potential K defined by (j3.1Uj) . that is: 

y^i^) = I ^j^) ''''''' II ^ S 

e/(a;), |a;| > 1 
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where / is a non-negative, continuous and radially symmetric function on {|x| > 1}. 
Moreover we assume / to be subcritical, that is it satisfies (jSHH): 

f{x) < K\x\-^-^, \x\>l, 

for some a,K > 0. 

We retain the notation of section 3.1, and in particular we recall the definition ()3.13p 
of eo- We have 

Theorem 4.3 (the operator —A — K on R^) For any e G (0, eo) the heat kernel of the 
operator —A — Ve satisfies 

K{t,x,y) < ct-^/^max{\x\~^,l}meix{\y\-^,l}, t>0, x,y eR^. (4.13) 



Proof. Let ip (x) be as in the proof of Theorem 13.41 cf ()3.15|) . It follows from ()3.14|) that 

r / I- \ (N-2)/N 

/^^ |V^| Vc^^ > c (/^Jt;p^/(^-^V'^/^^-')^r"'^^^''"'^rf^^ , (4.14) 

for all V e C](R^ \ {0}). Since xi2N-2)/iN-2) ^ ^2 ^ji^-^n y |o}) is a form 

core for and we conclude that K^{t, x, y) < ct~^/'^ whence, 

K{t,x,y)<ct-'"^,lj{x)ij{y). 

The required estimate on K(t,x,y) follows if we note that 

ci max{|a;|-(^-2)/^ 1} < tPix) < max{|x|-(^-2)/^ 1} , x G R^^. 

// 

It is well known that the estimates of the above theorems can be improved to yield 
Gaussian decay of the heat kernel. We have 

Proposition 4.4 Proposition \4 1\ as well as Theorems \4-^ and \4-^ can be improved by 
adding a factor cs exp{ — |x — yp/ ((4 + S)t)} to the corresponding right-hand sides. 

Proof. The proof is standard. One can use Davies's method of exponential perturbation 
[D] or Theorem 1.1 of lUj. Note that the argument is applied to the operator Hfp - not to 
H. ^ II 



5 Logarithmic refinements 

Our aim in this section is to obtain refined versions of the improved Hardy-Sobolev inequal- 
ities of Section 3. As an application, we prove heat kernel estimates for H = —A — V — Vi 
where V is one of the potentials studied in Section 4 (that is, V{x) = {{N — 2)/2)^|a;|~^ 
near zero) but Vi > is also critical. The criticality of Vi is meant in the sense that the 
following improved Hardy inequality holds 

/ ivwp - [ Vu^> [ Viu^, u G c^in), 

Jn Jn Jn 
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whereas this inequahty is no longer true if we replace Vi by (1 + e)Vi. Of course, Vi is of 
lower order with respect to \x\~^ (near x = 0) since ((A^ — 2)/2)^|x|~^ is already critical 
for the validity of the (simple) Hardy inequality. It is remarkable that the addition of the 
extra potential Vi does not affect the time dependence of the heat kernel estimates, but 
only affects the spatial singularity at the origin, which is increased by a logarithmic factor; 
see Theorems 15.31 and 15.41 

More precisely, recaUing that Xi{t) = (1 — logt)^^, let us introduce the functions 

X,+i(t) =Xi(X,(t)), fc = l,2,..., tG(0,l). (5.15) 

These are iterated logarithmic functions that vanish at an increasingly slow rate at t = 
and are equal to one at t = 1. In |FTj the following improved Hardy inequality was 
obtained for a bounded domain Q with D = supQ \x\ : 

1 /" ?/ ^ I T* I I I 

The potentials in the left-hand side of (j5.16j) are critical for each k, in the sense that (j5.16p 
is sharp: the term cannot be replaced by c^Xl^l for any e > 0, and the constant 1/4 
in the right-hand side is also optimal. In Theorem 15 . 31 and for bounded Q we obtain upper 
estimates on the heat kernel of the operator 

H--^- i^fj^ - I X? ^ ^ ^ - J^lX? ^ ^ ^ (5 AT) 

for /X > 1/4, as well as for the critical case fi = 1/4; for this we use results obtained in 
jFTj . For the critical case /i = 1/4 we also consider operators defined on R^, in analogy 
to the operator —A — K of Theorem 14. 3t for this we use Theorem 15.11 below, and the 
corresponding heat kernel estimate is given in Theorem 15.41 



5.1 Refined Hardy-Sobolev inequalities 

In this subsection we prove two theorems that are refined versions of Theorems 13.41 and 
EiH correspondingly. We recall definition ()5.15|) and set 

Y (Irr-W I ^k{\x\), \X\ < 1, 
^^^'"^'^ - \ 1, \X\ > 1, 

r n,(|x|) = Xfe(i/|x|), |x|>i, 

\ Yk{\x\)=Xk{l/\x\), \x\>0. 
We point out the differentiation rules for Xk{r) and Yk{r): 

^X^, = ^X,...X,.,Xr\ ^Y: = -'^Y,Y,...Y,^,Yr\ r = \xl (5.18) 
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valid for < r < 1 and r > 1 respectively, which are easily proved by induction. 

As in Theorem 13. 4t we assume that / is a non-negative, continuous and radially sym- 
metric function on B'^ satisfying 1)3.111) . that is, 

f{x) < K\x\-^-'', \x\>l, 

for some a,K > 0. For e > we also define 

\ e/(x), |x| > 1. 

We then have 

Theorem 5.1 Assume that k < N — 2 and define 

u€m(B'=) fu^dx 
Then ekfi > and for e G (0, ek,o) there holds 

iVnpda; - J^^ Vk,,u^dx > c (^J^^ . . . Xk+i)'^dx^^''~'^^'', (5.21) 



for all u e C~(R^) 

Remark. The constant e^^o is optimal in the sense that inequality ()5.21|) fails for e = e^^o- 
Also the exponent (2A^ ~ 2)/ (A^ — 2) in ()5.21|) is sharp in the sense that it cannot be 
replaced by a smaller exponent. The proof of these two facts is rather involved; see [FT] 
for similar arguments. We do not use these facts in the sequel. 

Proof. The proof follows closely that of Theorem 13. 4[ so we only give a sketch of it. The 
positivity of ekfi follows from Lemma EiH (i), yielding ek,o > K~^jiB{{,N — 2 + k) /2). Now 
let e G (0, Cfc^o) be fixed and let %Ij be the radially symmetric solution to the problem 



At/- + Vfc,,7/' = 0, |x|>l, 
normalized so that if) = 1 on = 1}. The function 



is then C^, radially symmetric and a direct computation which uses ()5.18p shows that 
A^ + Vk,eip = in R^. Exactly as in Theorem 13. 41 ip is positive, radially symmetric and has 
a positive limit as r — +oo. We then prove ()5.21|) in the case where u is radially symmetric, 
usmg once agam Proposition O The vahdity of ^1T\ for general u G C~(R^) follows 
from Lemma 13.31 / / 
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We finally prove a refined version of Theorem 13 .51 Let us fix a a non-negative, contin- 
uous and radially symmetric function g on B = {\x\ < 1}, such that 

g{x) < Klxl-^^", \x\<l, 

for some a,K > 0. Further for e > we define 

eg{x), \x\ < 1, 



+ |x|>i. 

We then have 
Theorem 5.2 Assume that k < N — 2 and define 

u<^m(B) j^gu^dx 
Then tkfl > and for e G (0, ek,o) there holds 

J^^ \Vu\^dx - J^^ Vk,,u^dx > c (1^^ . . . n+O^t^a;^^"'^^'^ , (5.23) 

for allue C^iYl^). 

Proof. We omit the proof, which is similar to that of Theorem I3.51 / / 

5.2 Refined heat kernel estimates 

In Theorems 14.21 and 14.31 we obtained heat kernel estimates for operators —A — V where 
V{x) = {{N — 2)/2)^|x|~^ near the origin. We shall now prove estimates for —A — V — Vi, 
with Vi also critical near the origin. In Theorems 15.31 and 15.41 we consider the cases ft 
bounded and Q = respectively. 
For k >1 and /i < 1/4 we define 

V,Ax) = ^ + ^,EXl..X! + ...XI xe n, (5.24) 

{Xi = Xi{\x\/D), D = supf^ |x| ) and consider the operator H = —A — V^f subject to 
Dirichlet boundary conditions on dQ. In |FT| Proposition 7.2] the Hardy-Sobolev inequal- 
ity 

/jlV^P - V,'^\')dx > c{J^ |nr^/(^-2)(Xi . ..X.+.f-^dxy'''"^^'', u e C^in), 

(5.25) 

was obtained. Let P be the largest solution of /?(1 — /?) = /i and define 

cPkA^) = \x\-^X['/' . . . Xl]I^Xl^ (5.26) 

Using (|5.18|) we verify that A(f)k^i3 + V^(j)k,/3 = and hence the change of variables u = 4>k,f3W 
yields 

/ (|Vnp-V;V)rfx= / \Vw\'<f>lpdx (5.27) 
for all w E C~(^]). We have 
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Theorem 5.3 Let Q be bounded, 1 < k < N — 2, and < /i < 1/4. The heat kernel of 
H = —A — satisfies the estimate 

K{t, X, y) < ct-^/^(l)kA^)(l)kAy), t>0, x,yen; (5.28) 

here Vj^ is given by \5.24\l and (pk,f3{x) by \5.2(^) . 

Proof. For the proof we distinguish two cases. 
(1) Case /i < 1/4. For w G C'^{^) we have 

/ \S/w\^(l)lgdx = [ \Vu\'^dx- [ V^u^dx 
Jn Jn Jn 

(bydnnD) > c(^J^{\Vu\^-Vjl!_y)dx^ 

(by > c(2 |^P^/(^-^)(^i • ..X,)'^dxy''-'^^'' 

{N-2)/N 



n 



This implies that K^^^{t,x,y) < ct'^/'^ and (jK^ follows. 

(2) Case fi = 1/4. By |FT1 Lemma 7.5] the following Sobolev inequality holds 



/ I VwHxI^-^Xf ^ . . . X^^dx > 
Jn 



2JV-2 



{N-2)/N 



> c[jjw\—^\x\-''Xi...XkXfJi-^'dxy , wec^in). 

This implies in particular 

/ \vw\'<plpdx>c{f \w\^^<pl^dxy''~'^^'', wec^in) 

and hence we have the uniform estimate K^^ ^(t, x, y) < ct'^/"^ as required. / / 

We finally have the following consequence of Theorem 15. H where we retain the notation 
of that theorem: 

Theorem 5.4 Let 1 < k < N — 2, and e G (0, ek,o) with ek,o given by ^5.2(^) . the heat 

kernel of the operator — A — Vk^t satisfies 

K{t, X, y) < ct'^^^ip{x)ilj{y), t > 0, x,y e R^; 
here 14, ^ is given by \5.1f^) and ip{x) by \5. 2^) . 

Proof. The result follows directly from Theorem 15.11 by means of ()4.6|) . / / 



22 



References 

[AE] Adimurthi, Esteban M. J. An Improved Hardy-Sobolev Inequality in W^'^ and its 
Application to Schrodinger Operator. Preprint 

[BV] Brezis H. and Vazquez J.-L. Blow-up solutions of some nonlinear elliptic problems. 
Rev. Mat. Univ. Comp. Madrid 10 (1997) 443-469. 

[CKN] Caffarelli L.A., Kohn R. V. and Nirenberg L. First order interpolation inequalities 
with weights. Compositio Math. 53 (1984) 259-275. 

[CM] Cabre X. and Martel Y. Existence versus instantaneous blowup for linear heat 
equations with singular potentials. C.R. Acad. Sci. Paris Ser. I Math. 329 (1999) 
973-978. 

[D] Davies E.B. Heat kernels and spectral theory. Cambridge University Press 1989. 

[DD] Davilla J. and Dupaigne L. Comparison Results for PDE's with a Singular Poten- 
tial. Preprint 2001. 

[DS] Davies E.B. and Simon B. norms of non-critical Schrodinger semigroups. J. 
Funct. Anal. 102 (1991) 95-115. 

[FT] Fihppas S. and Tertikas A. Optimizing Improved Hardy inequahties. J. Funct. Anal. 
192 (2002) 186-233. 

[G] Grigory'an A. Gaussian upper bounds for the heat kernel on arbitrary manifolds. 
J. Diff. Geom. 45 (1997) 33-52. 

[LN] Li Y. and Ni W.-M. On conformal scalar curvature equations in R^. Duke Math. 
J. 57 (1988) 895-924. 

[M] Maz'ya V. Sobolev spaces. Springer 1985. 

[MS] Milman P. and Semenov Yu. Heat kernel bounds and desingularizing weights. J. 
Funct. Anal., to appear. 

[T] Tertikas A. Critical phenomena in linear elliptic problems. J. Funct. Anal. 154 
(1998) 42-66. 

[VZ] Vazquez J.L. and Zuazua E. The Hardy inequality and the asymptotic behaviour 
of the heat equation with an inverse-square potential. J. Funct. Anal. 173 (2000) 
103-153. 

[WZ] Wang Zhi-Qiang and Zhu Meijun Hardy Inequalities with Boundary Terms. 
Preprint. 

[Z] Zhang Q.S. Global bounds of Schrodinger heat kernels with negative potentials. J. 
Funct. Anal. 182 (2001) 344-370. 



23 



